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ABSTRACT 

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel’dovich-selected 
galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sam- 
ple used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance 
end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We 
focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cos- 
mology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with 
a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives 
<78 = 0.851 ± 0.1 1 5 and w = — 1 . 14 ± 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year 
priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 
7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation 
obtained from numerical simulations and calibrated by X-ray observations, we find <78 = 0.821 ± 0.044 and 
w = -1.05 ±0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations 
plus type la supemoava which give erg = 0.802 ± 0.038 and w = -0.98 ± 0.053. A stacking analysis of the 
clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. 

These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and 
the cosmological parameters derived from them are broadly consistent with current models. 

Subject headings: cosmic microwave background - galaxies: dusters: general - cosmology: observations 
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1. INTRODUCTION 

Ever-improving observations suggest a concordant picture 
of our Universe. In this picture, generally called ACDM, 
"dark energy," the component responsible for the Universe’s 
accelerated expansion, is believed to be the energy of the vac- 
uum with a constant equation of state parameter, w, equal 
to —1 (e.g., Riess et al. 2009; Brown et al. 2009; Hicken et al. 
2009; Kessler et al. 2009; Percival et al. 2010; Komatsu et al. 
2010). ACDM has been measured via probes of the Uni- 
verse’s expansion rate such as type la supemovae, the pri- 
mary cosmic microwave background, and baryon acoustic 
oscillations. However, ACDM also makes concrete predic- 
tions about the Universe’s growth of structure. This growth 
rate describes how quickly dark matter halos form and evolve 
over cosmic time. A deviation from this predicted growth 
rate, particularly on linear scales, would signal a breakdown 
of ACDM (see e.g., Linder (2005); Bertschinger & Zukin 
(2008); Silvestri & Trodden (2009); Jain & Khoury (2010); 
Shapiro et al. (2010) and references therein). 

A handful of techniques have been available for measur- 
ing the growth of structure in the Universe. These largely 
consist of observing the weak and strong lensing of back- 
ground sources by intervening matter (e.g., Schrabback et al. 
2010), measuring distortions in redshift space with spectro- 
scopic surveys of galaxies (e.g., Simpson & Peacock 20 1 0), 
and quantifying the abundance of galaxy clusters as a func- 
tion of mass and redshift (e.g., Bahcall & Fan 1998)). This 
latter technique is one of the oldest and has been maturing 
with the advent of large area X-ray (e.g., Truemper 1990) and 
optical (e.g., Koester et al. 2007) surveys. 

Millimeter-wave surveys now possess the resolution and 
sensitivity to detect galaxy clusters. Detecting galaxy clus- 
ters via the Sunyaev-Zel’dovich (SZ) effect in large area 
millimeter-wave maps, as have become available through the 
Atacama Cosmology Telescope (ACT) (Swetzetal. 2010) 
and the South Pole Telescope (SPT) (Carlstrom et al. 2009), is 
a potentially powerful method. Cluster selection using the SZ 
effect (ZeTdovich & Sunyaev 1969; Sunyaev & ZeTdovich 
1970, 1972) is the technique whose selection function is least 
dependent on cluster redshift. This allows for a complete pic- 
ture of the evolution of clusters from their first formation to 
the present. 

Here we probe structure growth with a measurement of the 
abundance of massive galaxy clusters from observations made 
by the ACT project in 2008. We focus on the most massive 
SZ-selected clusters as this is the regime where high signal- 
to-noise measurements exist and we can best understand the 
cluster astrophysics. We also note that the clusters consid- 
ered in this work are rare and represent the tail of the mass 
distribution, which is sensitive to the background cosmol- 
ogy. This analysis uses the number of massive galaxy clus- 
ters to constrain, in particular, the normalization of the matter 
power spectrum, erg, and the dark energy equation-of state- 
parameter, w. 

This paper is structured as follows: Section 2 describes the 
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SZ effect and the ACT SZ cluster survey. Section 3 describes 
the 2008 ACT high-significance cluster sample. In Section 
4, we present our results, and in Section 5, we discuss their 
implications and conclude. 


2. BACKGROUND 

2.1. The Thermal Sanyaev-Zel 'dovich Effect 

The thermal SZ effect arises when primary cosmic mi- 
crowave background photons, on their path from the last 
scattering surface, encounter an intervening galaxy cluster. 
The hot ionized gas within the cluster inverse Compton 
scatters about 1% of the CMB photons, boosting their en- 
ergy and altering the intensity of the microwave background 
as a function of frequency at the location of the cluster 
(Sunyaev & ZeTdovich 1970, 1972). To first order, the effec- 
tive temperature shift (which is proportional to the intensity 
shift), at a frequency v. from the thermal SZ effect is given by 


AT 

Tcmb 


= /(*) 


ks&T 
m e c 2 


/ 


nJAl = f(x)y , 


( 1 ) 


where n e and T e are the number density and temperature of 
the electron distribution of the cluster gas, d / is the line-of- 
sight path length through the cluster, ay is the Thomson cross- 
section, kg is the Boltzmann constant, and 


fix) = xcoth(.r/2)-4, x = hv/(k^Tcmf- (2) 

Herey is the usual Compton y-parameter. Note that the full 
SZ effect contains relativistic corrections as in Nozawa et al. 
(1998), which we include in our simulations. For our sample, 
these corrections are 5 to 10%. We take Eq. 1, describing 
the first-order thermal SZ effect, as the definition of y we use 
throughout this work, and treat the relativistic corrections as 
an additional source of noise (see Section 3.3). 

At frequencies below 218 GHz, where the signal is null, 
AT is negative, and the cluster appears as a cold spot in CMB 
maps. Above the null, AT is positive, and the cluster ap- 
pears as a hot spot. Eq. 1 is also redshift independent, and the 
amplitude of the intensity shift is to first order proportional 
only to the thermal pressure of the cluster. This makes the 
SZ effect especially powerful for two reasons: the microwave 
background can trace all the clusters of a given thermal pres- 
sure that have formed between the last scattering surface and 
today in a redshift independent way, and the amplitude of this 
effect, being proportional to the thermal pressure, is closely 
related to the cluster mass. 


2.2. The ACT Sunyaev-ZeT dovich Cluster Survey 

The Atacama Cosmology Telescope (ACT) is a 6-meter off- 
axis telescope designed for arcminute-scale millimeter-wave 
observations (Swetz et al. 2010; Hincks et al. 2009). It is lo- 
cated on Cerro Toco in the Atacama Desert of Chile. One 
goal of this instrument is to measure the evolution of struc- 
ture in the Universe via the SZ effect. In the 2008 observ- 
ing season ACT surveyed 455 square degrees of sky in the 
southern hemisphere at 148 GHz. In this survey, galaxy clus- 
ters were detected from their SZ signal (see Marriage et al. 
(2010a) for details). A sample of 23 SZ-selected clusters 
was optically confirmed using multi-band optical imaging 
on 4-meter telescopes during the 2009B observing season 
(see Menanteau et al. (2010a) for details). Some of the low- 
redshift systems in this sample are previously known clus- 
ters for which spectroscopic redshifts are available. However, 
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roughly half are newly detected systems, and photometric red- 
shift estimates have been obtained from optical imaging. Here 
we make use of the subsample of these clusters with high- 
significance SZ detections (signal-to-noise ratio > 5, as de- 
fined in Section 3.2) to obtain cosmological parameter con- 
straints. 

3. THE 2008 ACT HIGH-SIGNIFICANCE CLUSTER 
CATALOG 

3.1. C MB Data 

Here we give a brief overview of the survey observations 
and the reduction of the raw data to maps. For a more com- 
plete introduction to the ACT instrument, observations, and 
data reduction pipeline, we refer the reader to Fowler et al. 
(2010) and Swetz et al. (2010). The 2008 observations in the 
southern hemisphere were carried out between mid-August 
and late December over a 9° wide ACT strip centered on 
a declination of -53° degrees and extending from approxi- 
mately 19* through 0* to 7*30 in right ascension. The 455 
square degrees used for this analysis consists of a 7°-wide 
strip centered at a declination of -52°30' and running from 
right ascension 00* 12” to 7* 1 0" ! . The resolution of the ACT 
instrument is about 1.4' at 148 GHz. Typical noise levels in 
the map are 30 pK per square areminute, rising to 50 /iK to- 
ward the map boundaries. Seven of the nine clusters consid- 
ered in this work fall in the central region of the map with 
lower noise levels. 

Rising and setting scans cross-link each point on the sky 
with adjacent points such that the data contain the information 
necessary to make a map recovering brightness fluctuations 
over a wide range of angular scales. In addition to survey ob- 
servations, ACT also executed regular observations of Uranus 
and Saturn during 2008 to provide beam profiles, pointing, 
and temperature calibration. Analysis of the beam profiles 
is discussed in Hincks et al. (2009). Absolute pointing is de- 
termined by comparing the positions of ACT-observed radio 
sources with the positions of these same sources detected in 
the AT20G survey (Murphy et al, 2010). The final tempera- 
ture calibration at 148 GHz is based on a recent analysis cross- 
correlating ACT and WMAP maps and is determined with an 
uncertainty of 2% (Hajian et al. 2010). The small residual 
calibration uncertainty translates into a small systematic un- 
certainty in y values of observed clusters, which is negligible 
compared to other uncertainties discussed in this analysis. 

To make maps, an iterative preconditioned conjugate gradi- 
ent solver is used to recover the maximum likelihood maps. 
This algorithm solves simultaneously for the millimeter sky 
as well as correlated noise (e.g., a common mode from at- 
mospheric emission). The map projection used is cylindrical 
equal area with a standard latitude of-53°30 r and square pix- 
els, 0.5' on a side. 

3.2. Cluster Detection Method 

In order to detect clusters in single-frequency millimeter- 
wave maps we construct a filter that is similar in morphology 
to the clusters we are trying to detect. We adopt a matched 
filter of the form 





]r(kQ [ 2 

P( k') 



r(k) 
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( 3 ) 


following Haehnelt & Tegmark (1996), Herranzetal. 
(2002a,b), and Melin et al. (2006). Here r(k) is the beam 


convolved cluster signal in Fourier space, and P( k) is the 
power spectrum of the noise, both astrophysical and instru- 
mental. The astrophysical noise sources for cluster detection 
include the primary CMB lensed by intervening structure, 
radio galaxies, dusty star-forming galaxies. Galactic dust, and 
the SZ background from unresolved clusters, groups, and the 
intergalactic medium. Since the power from the SZ signal is 
subdominant to these astrophysical sources (as evidenced by 
Luekeretal. (2010), Hall et al. (2010), Fowler etal. (2010), 
Dasetal. (2010), and Dunkleyetal. (2010)), we can to a 
good approximation model the pow'er spectrum of the total 
noise as the power spectrum of the data itself. In Eq. 3, the 
quantity in square brackets serves as a normalization factor 
to ensure an unbiased estimate of the cluster signal. When 
multi-frequency maps are available, this filter can be modified 
to incorporate the known spectral signature of the SZ signal. 

The template shape that we choose to match the cluster 
morphology is given by a two-dimensional Gaussian profile, 
which in Fourier space has the form 1 

AT{l) = Aexp[-8 2 (l+\)l/2]. (4) 

Here 6 = FWHM / \/81n2 where FWHM is the full width at 
half maximum, and A is a normalization factor that will be de- 
rived from simulations (see Section 3.3). We choose FWHM 
to be 2' as this is a typical cluster size in our maps. The analy- 
sis presented here of the cosmological parameters is nearly in- 
dependent of the particular profile chosen for the cluster tem- 
plate, as long as the template is smooth and well-matched to 
the cluster angular size. 

Before filtering our map to find clusters, we multiply the 
map, pixel-wise, by the square root of the number of obser- 
vations per pixel normalized by the observations per pixel in 
the deepest part of the map in order to establish uniform noise 
properties. We then detect point sources (radio and infrared 
galaxies) by a matched filter with the ACT beam as the tem- 
plate. Selecting all point sources with a signal-to-noise ratio 
greater than 4.0 in this filtered map, we mask them by re- 
placing all on-source pixels with signal-to-noise ratio greater 
than 4.0 with the average of the brightness in an annulus 4' 
away from the source center. We do this to avoid false de- 
tections due to the filter ringing around bright sources. See 
Marriage et al. (2010b) for details regarding point source de- 
tection. 

After masking out the brightest point sources, we filter the 
map to find clusters. Clusters are then detected within this 
filtered map with a simple peak detection algorithm along 
the lines of SExtractor (Bertin & Amouts 1996). An SZ y 
value for the brightest 0.5' pixel is measured for each clus- 
ter using this filtered map. This definition of v is different 
from the integrated Y, which is specifically the Compton-y 
parameter integrated over the face of the cluster and given by 
Y = J yd fl. The integration for this Y value is performed over 
a radius tied to the size of the cluster, and a Y defined this 
way would be a preferable quantity to use, having lower scat- 
ter with mass in theory (e.g., da Silva et al. 2004; Motl et al. 
2005; Nagai 2006; Reid & Spergel 2006; Bonaldi et al. 2007). 
However, given single-frequency millimeter-wave maps, the 
size of each cluster cannot always be robustly determined. An 
alternative quantity to measure is a "central y value," gener- 
ally referred to as yo, which essentially describes the normal- 
ization of the specific template shape used to find the cluster. 

1 We use the flat space approximation, I = 2trk. 
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TABLE 1 

ACT Cluster Catalog for High-Significance Clusters from the 2008 Observing 

Season 


ACT Descriptor 

R.A. 

dec). 

y7cMB(/iK) 

Redshift 

Other Name 

ACT-CI. J0645-5413 

06:45:30 

-54:13:39 

340 ± 60 

0.1 67 a 

Abell 3404 

ACT-CL J0638-5358 

06:38:46 

-53:58:45 

540 ± 60 

0.222 a 

Abell S0592 

ACT-CL J0658-5557 

06:58:30 

-55:57:04 

560 ± 60 

0.296 b 

1 ES0657-558(Bullet) 

ACT-CL J0245-5302 

02:45:33 

-53:02:04 

475 ± 60 

0.300 c 

Abell S0295 

ACT-CL J0330-5227 

03:30:54 

-52:28:04 

380 ± 60 

0.440 d 

Abell 3128(NE) 

ACT-CL J0438-5419 

04:38:19 

-54:19:05 

420 ± 60 

0.54±0.05 e 

New 

ACT-CL J06 16-5227 

06:16:36 

-52:27:35 

360 ± 60 

0.71 ±0.10 c 

New 

ACT-CL J0102-4915 

01:02:53 

-49:15:19 

490 ± 60 

0.75 ±0.04 e 

New 

ACT-CL J054 6-5345 

05:46:37 

-53:45:32 

310 ±60 

1.066' 

SPT-CL 0547-5345 


1 /iK given for the brightest 0.5 / pixel of each duster 
a spectroscopic-z from de Grandi et al. (1999) 
b spectroscopic-z from Tucker et al. (1998) 
c spectroscopic-z from Edge etal. (1994) 
d spectroscopic-z from Werner et al. (2007) 
e photometric-z from Menanteau et al. (2010a) 
f spectroscopic-z from Infante et al. (2010); Brodwin et al. (2010) 


This quantity is not ideal for a cosmological analysis as it is 
intimately tied to the profile shape whereby the yo value is 
determined. These values also exhibit a larger scatter with 
cluster mass than an integrated y quantity (e.g., Motl et al. 
2005). For all the clusters considered in this analysis, we fix 
an aperture size, given by our pixel size of 0.5', and mea- 
sure Compton-y values within this fixed aperture. We do not 
consider larger aperture sizes here because we wish to limit 
template shape and redshift dependence. 

We find that selecting clusters with a ,v7cmb value 2 greater 
than 300 //K corresponds to a subsample of clusters with 
a signal-to-noise ratio greater than 5. Here signal-to-noise 
ratio is defined as the signal of the brightest cluster pixel 
in the filtered map divided by the square root of the noise 
variance in the filtered map. This subsample corresponds 
to the subsample of clusters with signal-to-noise ratio >5.9 
in Marriage et al. (2010a). 3 The y values derived with the 
method used in this paper are not directly comparable to those 
in Marriage et al. (2010a), in which an optimal filter sized is 
searched for. However, the methods are independently com- 
pared to simulations. While we detect clusters down to a 
signal-to-noise ratio of about 3 as defined in Marriage et al. 
(2010a), we use only this higher-significance subsamplc in 
this work. This subsample is given in Table 1 . 

3 .3 . Simulations and SZ Signal Recovery 

To determine the expected scatter in our recovered vTcmb 
values, we perform the same detection procedure discussed 
above on simulated maps. Hereafter, the simulations we re- 
fer to are those discussed in Sehgal etal. (2010), which in- 

2 Note that y is a dimensionless parameter. We multiply it by 7 cmb = 
2.726 X 1 0° /(K to give an indication of the expected temperature decrements. 
For the frequency dependence, /(a) w -1 in Eq. 1 at 148 GHz. 

3 In Marriage et al. (2010a), a different detection method is used that varies 
the angular scale of the filter to match clusters of different sizes, and assigns 
a signal-to-noise ratio based on the scale that gives the highest value. The one 
duster that has a signal-to-noise ratio > 5,9 in that work that is not included 
here is ACT-CL J0235-5 121. In Marriage et al. (2010a), it was found to have 
a high signal-to-noise ratio using a template scale of 4.0'. Although there is 
no doubt that this is a massive cluster (Menanteau et al. 2010a), its relatively 
high redshift, z = 0.43 ± 0.07 argues for a compact size, suggesting that CMB 
contamination could be boosting the clusters ’s signal-to-noise ratio on a 4.0' 
scale, as discussed in Marriage et al. (2010a). This cluster is not found with 
signal-to-noise ratio > 5 using the 2' FWHM Gaussian template described 
above. 


elude the SZ signal, lensed primary cosmic microwave back- 
ground, Galactic dust, and radio and infrared sources cor- 
related with SZ clusters as suggested by observations. The 
large-scale structure in this simulation was carried out using 
a tree-particle-mesh code (Bode et al. 2000; Bode & Ostriker 
2003), with a simulation volume of 1000 /T*Mpc on a side 
containing 1 024 3 particles. The cosmology adopted is con- 
sistent with the WMAP 5-year results (Komatsu et al. 2009) 
though the details of the cluster properties are relatively in- 
sensitive to the background cosmology. The mass distribution 
covering one octant of the full sky was saved, and halos with a 
friends-of-friends mass above 1 x 10 I3 M Q and with a redshift 
below z = 3 are identified. The thermal SZ signal is derived by 
adding to the N-body halos a gas prescription that assumes a 
polytropic equation of state and hydrostatic equilibrium. This 
model, which is described in more detail in Bode et al. (2009), 
adjusts four free parameters (star-formation rate, nonthermal 
pressure support, dynamical energy transfer, and feedback 
from active galactic nuclei) which are calibrated against X- 
ray gas fractions as a function of temperature from the sample 
of Sun et al. (2009) and Vikhlinin et al. (2006). The pressure 
profiles of the massive, low-redshift clusters in this simula- 
tion agree well with the best-fit profile of Amaud et al. (2009) 
based on X-ray observations of high-mass, low-redshift sys- 
tems (Trac et al. 2010). We will see in Section 4.5 that the 
stacked SZ signal of the clusters in Table 1 is also consistent 
with the stacked thermal SZ signal of the massive clusters in 
this simulation. The kinetic SZ in this simulation is calculated 
from the line-of-sight momentum of the particles. We also 
include the relativistic corrections to the SZ signal as given 
in Nozawa et al. (1998). We convolve these simulations with 
the ACT beam and run them through the same map-making 
process discussed in Section 3.1, including simulated atmo- 
spheric emission and realistic instrumental noise. 

From these simulations, we cut out six different patches of 
455 square degrees to mimic the sky coverage in this analysis. 
These six sky patches give us about 40 clusters that would cor- 
respond to the high-significance cluster sample given in Table 
1 . Using these simulations, we apply the same cluster detec- 
tion procedure as discussed in Section 3.2, and recover v7cmb 
values for the detected clusters. These recovered vTcmb val- 
ues are compared to the true j7cmb values taken from the 
first-order thermal SZ maps alone, prior to any instrumental 
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Fig. 1 . — The truey value (as defined in Eq. 1) for the brightest 0.5' pixel 
of each cluster versus the recovered y value from simulations using the detec- 
tion method outlined in Section 3.2. The dimensionless y values have been 
multiplied by 7cmb = 2.726 X 10 6 /uK to give an indication of the expected 
temperature decrements at 148 GHz. The root-mean-square scatter shown 
here is 60 fj, K. 

or atmospheric modifications. The comparison between these 
true and recovered j;7cmb values is shown in Figure 1. We 
set the normalization factor, A, in Eq. 4 such that the mean 
bias between true and recovered v7cmb values is zero. The 
root-mean-square scatter in the recovered yTcm values is 60 
pK. The scatter here is dominated by the instrumental and 
atmospheric noise sources in the map and not by the tech- 
nique itself. This same normalization is used for filtering the 
data map, and the same scatter is assumed. Figure 1 also 
shows that any bias due to boosting a cluster with an intrinsic 
vTcmb < 300 fiK to a value above 300 //K is far below the 
scatter. Vanderlinde et al. (2010) found that this signal boost- 
ing effect is at most 4% for clusters detected with a signal-to- 
noise ratio greater than 5. 


Fig. 2. — The redshift distribution of the 9 high-significance clusters listed 
in Table 1 compared to expectations from the mass function of Tinker et al. 
(2008) assuming all clusters above a given mass threshold have been detected. 
Here we assume the total area observed is 455 square degrees with uniform 
coverage, and the WMAP7+BAO+SN best-fit cosmology for a wCDM model 
(Komatsu et al. 2010) for the mass function. This figure suggests an effective 
mass threshold of our sample of ss 10.4 x 10 14 Mq (dashed line). For illus- 
trative purposes we show via dotted lines the expected redshift distributions 
for mass thresholds larger/smaller by 0.6 x IO'LWq. 

terns we provide photometric redshifts based on the NTT and 
SOAR imaging. Table 1 gives the mean photometric redshift 
for these clusters obtained by iteratively selecting galaxies 
photometrically classified as E or E/SOs within a projected 
radius of 500 kpc and redshift interval |Az| = 0.05. This 
was done to obtain a local color-magnitude relation for each 
color combination using a 3n median sigma-clipping algo- 
rithm. The uncertainties on the photometric redshifts of the 
three new clusters in Table 1 come from the weighted rms 
of the individual galaxies chosen as members. Details of the 
photo-z algorithm are given in Menanteau et al. (2010b). 


3 .4. Optical Identification and Cluster Redshifts 

The sample of SZ-selected clusters obtained from the data 
via the method above was followed up with optical obser- 
vations to verify the millimeter-wave cluster identifications 
and determine cluster redshifts. Here we provide a summary 
of the observing strategy and redshift determinations of our 
cluster sample, and we refer the reader to Menanteau et al. 
(2010a) for a detailed description. 

SZ cluster candidates were observed during the 2009B ob- 
serving season with optical imaging on the 4-meter SOAR and 
NTT telescopes to search for a brightest cluster galaxy and an 
accompanying red sequence of cluster members. While some 
of the clusters in the SZ-selected sample correspond to previ- 
ously known systems at low redshift (z< 0.3), some represent 
new systems, previously undetected at other wavelengths. 

Photometric redshifts and their probability distributions, 
p(z), were computed for each object from their dust-corrected 
gri isophotal magnitude using the BPZ code (Benitez 2000). 
Six of the clusters in our sample have spectroscopic redshift 
information available (see Table 1 for references), as they 
were previously known systems. For the remaining three sys- 


4. RESULTS 

4.1. SZ Selection Function 

We determine the selection function of our cluster subsam- 
ple through both optical observations and simulations. To 
investigate the effective mass threshold of our cluster sam- 
ple, we plot cumulative clusters as a function of redshift 
and compare that to expectations from the mass function of 
Tinker et al. (2008) assuming the best-fit cosmology from 
WMAP7+BAO+SN with a wCDM model (Komatsu et al. 
2010). We find that this sample is consistent with a mean 
mass 4 threshold of 10.4 x 10 14 Mq for M200 as shown by the 
dashed line in Figure 2. For illustrative purposes, we show 
by dotted lines the expected redshift distributions for mass 
thresholds larger and smaller by 0.6 x 10 14 M©, which enclose 
our subsample distribution for z > 0.25. This figure also as- 
sumes uniform coverage over an area of 455 square degrees. 
Note that Figure 2 is presented to give a qualitative under- 

4 Note that throughout this text cluster masses are defined in terms of ATjoo. 
which is the mass within Ri oo, the radius within which the mean cluster den- 
sity is 200 times the average density at the cluster redshift. 
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Fig. 3. — Completeness versus true jFcmb from simulations. For an ob- 
served >'7cmb threshold of 300 /rK we expect our sample to be about 85% 
complete for true y7cMB values above 300 /rK. 


Fig. 4. — Purity versus recovered v7cmb from simulations. For an observed 
v7cmb threshold of 300 pK we expect our sample to be about 90% pure. Our 
actual sample is 1 00% pure, as each candidate cluster has been confirmed by 
optical observations (Menanteau et a], 2010a). 


standing of our sample, but is not directly used in the cosmo- 
logical analysis. 

We also use simulations to characterize our sample. The 
same cluster detection procedure discussed in Section 3.2 is 
applied to the simulated millimeter-wave maps discussed in 
Section 3.3, including instrumental and atmospheric noise 
sources. We find for an observed >’Tcmb threshold of 300 //K, 
we expect from simulations to be about 85% complete above 
a true j'7cmb value greater than 300 /iK as shown in Figure 
3. This completeness is calculated as the total number of ob- 
served clusters with a recovered y value above the threshold 
versus the total number of expected clusters with a truey value 
above the threshold. This takes into account the scatter of ob- 
served y values across the threshold. We take 300 //K as our 
threshold v7cmb value and only consider clusters with mea- 
sured vTcmb values larger than this. We expect about 90% 
purity for cluster detections with an observed jTcmb value 
greater than 300 /«K. Figure 4 illustrates how this false detec- 
tion rate is expected to vary as a function of observed i'7'cmb 
threshold. From the optical observations discussed above, we 
find that all clusters identified in the millimeter- wave maps 
with >’7 cmb values greater than 300 //K were verified as clus- 
ters in the optical. Thus this sample is 100% pure with no 
false detections (see Menanteau et al. (20 1 0a)). Note that Fig- 
ures 3 and 4 are representative, coming from a relatively small 
sample of simulated high-significance clusters as discussed in 
Section 3.3. They are presented to give a qualitative under- 
standing of the sample and to understand above what >’7cmb 
threshold our sample is roughly complete. Beyond this, they 
do not enter in the analysis of cosmological parameters. 


4.2. Scaling Relation Between SZ Signal and Mass 

We assume for the relation between the SZ signal and mass 
the general parameterized form 5 


y* me =A 



(5) 


5 Note that the form of this relation is analogous to that in 
Vanderlinde et al. (2010). 


where M 0 = 5 x IO^Mq/T 1 and zq = 0.5. We also assume 
that cluster y values are randomly distributed around this 
relation with a lognormal scatter S. By lognormal scatter, 
we mean that the scatter in the relation ln(y ,rae ) = \n(A) + 
B\n(M mc / M Q )+C\n((\ +z)/( 1 +zo)) has a Gaussian distribu- 
tion with a mean of zero and a standard deviation equal to 5. 
This scaling relation relates the observable quantity for each 
cluster, SZ signal, to the quantity most directly relevant for 
cosmology, the cluster mass. This form is chosen as the SZ 
signal is expected to have a power-law dependence on mass. 
Some dependence of our recovered SZ signal on redshift is al- 
lowed with an additional parameter C. Note that the B and C 
parameters are independent of Mq and zq, whose values only 
affect A. We also note that the uncertainty in the relation be- 
tween SZ signal and mass is dominated by astrophysical pro- 
cesses and is only minimally dependent on cosmological pa- 
rameters given WMAP7 priors. 

We use the simulations discussed in Section 3.3 to deter- 
mine fiducial model values for A, B, and C. In particular, we 
use only the simulated maps containing the first-order thermal 
SZ component of the simulated clusters, without altering the 
maps by adding any noise or convolving with the ACT beam. 
Using they value of the brightest 0.5' pixel for each cluster in 
the simulated thermal SZ maps, we solve for the best-fit val- 
ues of the three scaling relation parameters in Eq, 5, as well 
as the scatter in this relation, with a linear least squares fit. 
Figure 5 shows the simulated clusters along with the best-fit 
line, and Table 2 gives the corresponding best-fit values and 
errors of the fiducial scaling relation parameters in addition to 
the scatter S. The best-fit values for the B and C parameters 
are close to what we would expect from self-similar scaling 
relations as discussed in Appendix A. 

There are a number of astrophysical mechanisms that could 
cause the observed relation between SZ signal and mass to 
differ from the fiducial relation. One is contamination of 
SZ decrements by radio or infrared galaxies at 148 GHz. 
Regarding radio galaxies, observations suggest that these 
galaxies show some preference for residing in galaxy clus- 
ters (Coble etal. 2007; Lin & Mohr 2007; Linetal. 2009; 
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3 5 10 15 20 30 

Moo (10 14 M s h~ l ) 

FlO. 5. — Relation between true y and the cluster mass from simulations, 
including clusters with A/joo > 3 X IO^A/qA' 1 and z > 0.15. The best-fit 
scaling relation parameters in Eq. 5 are found with a least-squares fit, and 
the resulting best-fit surface is plotted in two-dimensions as the solid line 
above. The lognormal scatter between the truer values and the best-fit scaling 
relation is 26%. 

Mandelbaum et al. 2009). However, using a model of ra- 
dio galaxies that describes their correlation with halos, the 
amount of contamination expected from radio galaxies was 
found to be negligible for (Sehgal et al. 2010). For redshifts 
< 1, star formation, which is responsible for infrared galaxy 
emission, is expected to be quenched in high-density environ- 
ments. At low redshifts (z ~ 0.06) the fraction of all galaxies 
that are star fonning galaxies is ~ 16% in clusters (Bai et al. 
2010). While this percentage is expected to increase at higher 
redshifts, given that the total infrared background at 150 GHz 
is roughly 30 // K (Fixsen et al. 1998), it is unlikely that in- 
frared galaxy contamination could be significant for clusters 
with y7cMB > 300 fi K. Lima et al. (2010) have also shown 
that the lensing of infrared galaxies by massive clusters should 
not introduce a significant bias in the measured SZ signals. 

Another way for the observed SZ signal to be lower than 
the fiducial model is if clusters have a significant amount of 
nonthermal pressure. This pressure would not be observed 
as part of the SZ signal, however, it would play an important 
role in counteracting the gravitational pressure from the clus- 
ter mass. Such nonthermal pressure can take the form of small 
scale turbulence, bulk flows, or cosmic rays. Simulations 
and observations suggest contributions to the total pressure 
from cosmic rays to be about 5 - 10% (Jubelgas et al. 2008; 
Pfrommer & Enfilin 2004) and from turbulent pressure to be 
between 5-20% (Lauetal. 2009; Meneghetd et al. 2010; 
Bums et al. 2010), with only the latter work suggesting levels 
as high as 20% and that largely at the cluster outskirts. These 
processes have a much larger impact on lower mass clus- 
ters and groups where the gravitational potential is not strong 
enough to tightly bind the cluster gas (e.g., Battaglia et al. 
2010; Shaw et al. 2010; Trac et al. 2010). However, for the 
massive systems considered here, this again is not expected to 
be a significant issue. One astrophysical process that can have 
a significant affect on the cluster y values is major mergers. 
We certainly have at least one in our sample (Bullet cluster), 
but note the extreme rarity of such objects in general. 


4.3. Cluster Likelihood Function 

In order to constrain cosmological parameters with our 
cluster sample, we construct a likelihood function specific for 
clusters, and we map out the posterior distribution to find 
marginalized distributions for each parameter. We follow 
Cash (1979) who derived the likelihood function in the case 
of Poisson statistics giving 

Ni, 

In £ = InP )({«/}|{A,}) = ]T(n,lnA,- A,). (6) 

r=! 

P r is the probability of measuring {«;} given modeled counts 
{A,}. Here A* is the total number of observed bins in SZ sig- 
nal - redshift space, and A, is the modeled number of clusters 
in the /th bin. We also take the bin sizes to be small enough 
so that no more than one observed cluster is in each bin. The 
modeled cluster count. A,, is a function of the SZ signal and 
redshift of the given bin (which we call y obs and z obs ) as well 
as the set of cosmological parameters, {cj}. The modeled 
count is also a function of the parameters of the SZ signal - 
mass scaling relation ( A,B,C,S ) given in Eq. 5, since it is the 
abundance of clusters as a function of mass that is tied to cos- 
mology via the mass function. For this work we use the mass 
function given in Tinker et al. (2008). A derivation of the full 
cluster likelihood function used in this analysis can be found 
in Appendix B. This likelihood is given by Eq. B8 and is a 
function of the parameters {c ; } and d, TEC, S'. 

We assume normal errors of 2.2 x 1 0~ 5 on>> obs (correspond- 
ing to an error \nyT CM u of 60 fiK) and 0. 1 on z obs . We take 0. 1 
as the redshift uncertainty for convenience even though six of 
our clusters have spectroscopic redshifts. However, the red- 
shift error does not dominate the uncertainty of our results. 
We also assume Gaussian priors on A,B,C, and S centered 
around the fiducial values given in Table 2, with conservative 
la uncertainties of 35%. 20%, 50%, and 20% respectively of 
the fiducial values. These priors were determined by find- 
ing the relation between SZ signal and mass from simulated 
thermal SZ maps with varying gas models. In particular, we 
use two simulated thermal SZ maps analogous to those dis- 
cussed in Section 4.2, with the gas physics models in these 
maps based on the adiabatic and the nonthermal20 models de- 
scribed in Trac et al. (2010). The adiabatic model assumes no 
feedback, star-formation, or other nonthermal processes that 
could lower the SZ signal as a function of mass. The nonther- 
mal20 model assumes more star-formation than the fiducial 
model and 20% nonthermal pressure support for all clusters 
at all radii, which is a larger amount of nonthermal pressure 
than generally suggested by X-ray observations and hydro- 
dynamic simulations (e.g., Lau et al. 2009; Meneghetti et al. 
2010; Bums et al. 2010). These two models span the range 
of plausible gas models for massive dusters given current ob- 
servations, and the 1 a priors on the scaling relation parame- 
ters given above are generous given the range in parameters 
spanned by these models. 

4.4. Parameter Constraints 

The likelihood function described above was made into a 
standalone code module which was then interfaced with the 
Markov chain software package CosmoMC (Lewis & Bridle 
2002). Using CosmoMC, we run full chains for the WMAP7 
data alone (Larson et al. 2010) and for the WMAP7 data plus 
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FIG. 6. — Likelihood contour plots of ir versus a$ showing ler and 2a marginalized contours. Left-. Blue contours are for WMAP7 alone, and red contours are 
for WMAP7 plus ACT SZ detected clusters, fixing the mass-observable relation to the fiducial relation given in Section 4.2. Right: Contours are the same as in 
the left panel, except that the uncertainty in the mass-observable relation has been marginalized over within priors discussed in Section 4.3. 


TABLE 2 

Best-fit Scaling Relation Parameters 


Model and Data Set 

A 

B 

C 

S 

Simulation Fiducial Values 
wCDM WMAP7 + ACT Clusters 

(5.67 ± 0.05) X 1(T 5 
(8.77 ± 3.77) x 10" 5 

1.05 + 0.03 
1.75 + 0.28 

1.29 + 0.05 
0.97 + 0.68 

0.26 

0.27 + 0.13 


TABLE 3 

Cosmological Parameter Constraints for cr s and w 


Model and Data Set 


w 

wCDM WMAP7+BAO+SN 
wCDM WMAP7 

wCDM WMAP7 + ACT Clusters (fiducial scaling relation) 
wCDM WMAP7 + ACT Clusters (marginalized over scaling relation) 

0.802 + 0.038 
0.835 + 0.139 
0.821+0.044 
0.851+0.115 

-0.98 + 0.053 
-1.11+0.40 
-1.05 + 0.20 
-1.14 + 0.35 


our ACT cluster subsample. We assume a wCDM cosmo- 
logical model which allows w> to be a constant not equal to 
-1, assumes spatial flatness, and which has as free parame- 
ters: flbh 2 , £l c h 2 , 9 *, t, w, n s , ln[10 io zf,], and Asz as defined 
in Larson etal. (2010). The parameter ag is derived from 
the first seven of these parameters and is kept untied to As z 
as the link between the two is in part what we are investi- 
gating. We run the WMAP7 plus ACT clusters chain under 
two cases: one where the values of A,B.C , and 5 are fixed 
to the fiducial values given in Section 4.2 and listed in Table 
2, and one where A.B.C. and S are allowed to vary within 
the conservative priors given in Section 4.3. For the latter 
case, we add these four new parameters to the CosmoMC 
code. At each step of the chain, CosmoMC calls the soft- 
ware package CAMB 6 to generate both the microwave back- 
ground power spectrum and matter power spectrum as a func- 
tion of the input cosmology, and then the natural logarithms 
of both the WMAP and cluster likelihoods are added. We 
determine the posterior probability density function through 
the Markov chain process and use a simple R-l statistic 
(Gelman & Rubin 1992) of R-l < 0.01 to check for conver- 

6 www.camb.info 


gence of the chains. 

The best-fit marginalized ler and 2a contours, obtained 
from this process, are shown in Figure 6 for w and crs. 
The blue contours show the constraints for WMAP7 alone, 
while the red contours show the constraints from the union 
of WMAP7 plus our ACT cluster subsample. The left panel 
shows the best-fit contours with the SZ signal - mass scaling 
relation fixed to the fiducial relation obtained from the sim- 
ulations. The right panels show the constraints allowing the 
four parameters of the scaling relation to vary. Table 3 lists the 
best-fit parameter values for a% and w with their la marginal- 
ized uncertainties. Table 2 lists the best-fit scaling relation 
values and ler uncertainties as well as the fiducial values ob- 
tained from simulations as discussed in Section 4.2 for com- 
parison. We note that for the remaining seven parameters fit 
in the analyses combining WMAP7 plus ACT clusters (0*/t 2 , 
Q c h 2 , t, n s , ln[]0 l0 zf s ], and Asz X we find best-fit values 
consistent with the best-fit values from WMAP7 alone with a 
modest improvement in the marginalized errors. 

4.5. Slacked SZ Signal 

We also perform a stacking analysis of the nine clusters 
listed in Table I to measure average cluster SZ profiles, which 
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Fig. 7. — Average profile from stacking the 9 clusters presented in Table 
1 (black solid line) compared with the average profile stacking 40 clusters 
with unsmoothed j7cmb values greater than 300 juK in simulated thermal SZ 
maps convolved with the ACT beam (dashed blue line). The dashed blue line 
shows the average profile for clusters simulated with the fiducial SZ model, 
while the dotted red (bottom) and dot-dashed green (top) lines show the same 
assuming the adiabatic and nonthermal20 SZ models respectively which are 
discussed in Section 4.3, For the profiles of the 9 clusters in the data, we 
removed a mean background level from the profile of each cluster. Error bars 
for the simulated clusters have been offset by 0. 1' for clarity, and are smaller 
than those from the data. Error bars for the adiabatic and nonthermal20 mod- 
els are not shown, but are of similar size as for the fiducial model. 

can be compared with simulations. We stack the 9 clusters in 
the data map prior to any filtering, after subtracting a mean 
background level for each cluster profile using an annulus 15' 
from the center of each cluster and 0.5' wide. The stacked av- 
erage profile is given by the solid black line in Figure 7. The 
same procedure is preformed on all simulated clusters with 
unsmoothed yTcMB values greater than 300 //K in simulated 
thermal SZ maps. There are 40 of these simulated clusters 
in total over 6 different 455 square degree maps spanning the 
same redshift range as the data. These simulated clusters are 
stacked in thermal SZ maps convolved with the ACT beam 
to mimic the data, and their average profile is given by the 
dashed blue line in Figure 7. Error bars represent the standard 
deviation of the mean in each radial bin. The blue dashed line 
represents the stacked profiles of simulated clusters assuming 
the fiducial SZ model. The red dotted and green dot-dashed 
lines show the stacked profiles of simulated clusters assum- 
ing the adiabatic and nonthermaOO SZ models discussed in 
Section 4.3. The error bars have not been included for the lat- 
ter two models in Figure 7, but they are of similar size as for 
the fiducial model. We find good agreement in the average 
profiles of the clusters in the data and simulated with the fidu- 
cial model as shown in Figure 7, which suggests that there is 
no significant misestimate of the SZ signal for these massive 
systems. 

5. DISCUSSION 

From Table 3 we see overall agreement between cr$ and w 
as measured with only WMAP7 and as measured with the 
high-significance ACT cluster sample plus WMAP7. We find 
<78 = 0.821 ±0.044 and w = -1 .05 ± 0.20 if we assume the 


fiducial scaling relation, a decrease in the uncertainties on 
these parameters by roughly a factor of three and two re- 
spectively as compared to WMAP7 alone. This indicates 
the potential statistical power associated with cluster mea- 
surements. Marginalizing over the uncertainty in this scaling 
relation, we find erg = 0.851 ± 0.1 15 and w = -1.14 ± 0.35, 
an uncertainty comparable to that of WMAP7 alone. We 
also see consistency when comparing these constraints to the 
best-fit constraints from WMAP7 plus baryon acoustic oscil- 
lations plus type la supemovae, which give <78 = 0.802 ±0.038 
and w = -0.980 ± 0.053 for a wCDM cosmological model 
(Komatsu et al . 2010). As the latter are all expansion rate 
probes, this suggests agreement between expansion rate and 
growth of structure measures. Both also show w is consistent 
with -1 , giving further support to dark energy being an energy 
of the vacuum. 

These results are also consistent with analyses from X-ray 
cluster samples giving <r 8 (!n m /0.25) 0 - 47 = 0.8 1 3 ±0.013 (stat) 
±0.024 (sys) and w = -1 . 14 ± 0.21 (Vikhlinin et al. 2009), 
n m = 0.23 ±0.04, <7 S = 0.82 ±0.05, and w = -1.01 ± 0.20 
for a wCDM model (Mantz et al. 2010), and fl m = 0.30^'q2, 
<t 8 = 0.85^ q 2 f rom WMAP5 plus X-ray clusters (Henry et al. 
2009) . We also find consistency with optical samples yield- 
ing <7 8 (Q, n /0.25) 0 ' 41 = 0.832 ± 0.033 for a flat ACDM model 
(Rozo et al. 2010). Vanderlinde et al. (2010) find <7g =0.804 ± 
0.092 and w = -1 .049 ± 0.291 for a wCDM model using SZ 
clusters detected by SPT plus WMAP7. 

This analysis also suggests consistency between the fiducial 
model of cluster astrophysics used here to describe massive 
clusters and the data. Table 3 shows agreement in best-fit cos- 
mological parameters between growth rate and expansion rate 
probes when we hold fixed our fiducial relation between SZ 
signal and mass. When we allow the scaling relation parame- 
ters to be free, we find best-fit values that are broadly consis- 
tent with those of our fiducial relation. We note that while the 
1 <7 range of the B parameter is higher than the fiducial value, 
the fiducial value is enclosed by the 2a range of 1 .751 q 4 . The 
higher value of the B parameter may indicate some curvature 
in the true scaling relation away from the fiducial model at 
the high-mass end. This may also be suggested by Figure 5 
where the simulated clusters seem to prefer higher y values 
than the fiducial relation would suggest for the most massive 
systems. The agreement between cosmological parameters 
from expansion rate and growth of structure probes when fix- 
ing the SZ signal - mass scaling relation to the fiducial model 
and the broad agreement between fiducial and best-fit scal- 
ing relation parameters when the latter are allowed to be free, 
suggest our data is broadly consistent with expectations for 
the SZ signal of massive clusters. This is also suggested by 
comparing the stacked SZ detected clusters in the data with 
simulations as shown in Figure 7. 

We would expect the above to be the case as massive clus- 
ters have been studied far better than lower mass clusters with 
a variety of multi-wavelength observations. In addition, a 
number of astrophysical processes that are not perfectly un- 
derstood, such as nonthermal processes and point source con- 
tamination, affect the gas physics of lower mass clusters much 
more than that of the most massive systems. In general, these 
processes tend to suppress the SZ power spectrum over that 
of a straightforward extrapolation based on the most massive 
systems. This is an important effect as lower mass systems 
(< 10 14 M Q ) contribute as much to the SZ power spectrum at 
/ ~ 3000 as systems at higher mass (Komatsu & Seljak 2002; 
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Trac et al. 2010). The power spectrum near / = 3000 has been 
recently measured by and discussed in Luekeretal. (2010), 
Das et al. (2010) and Dunkley et al. (2010). 

There are a number of ways the cosmological constraints 
presented here could be further improved. Clearly the largest 
uncertainty is the relation between SZ signal and mass, and 
further X-ray observations of massive clusters, particularly at 
higher redshifts where X-ray observations have been limited, 
would help to calibrate this relation. Further targeted observa- 
tions of massive clusters at millimeter-wave frequencies with 
enough resolution and sensitivity to identify point sources 
would offer a better handle on contamination levels. In addi- 
tion, an analysis using multiple frequency bands, which would 
employ the spectral information of the SZ signal, may be 
helpful in determining cluster sizes and measuring integrated 
Ls. This could help reduce the scatter in the relation between 
SZ signal and mass. Spectroscopic redshifts of all the clusters 
in a given SZ sample would also help to reduce uncertainty 
on the cosmological parameters. In addition, millimeter- wave 
maps with lower instrument noise, would greatly reduce the 
scatter between the recovered and true SZ signal. Such maps 
are expected with ACTpol (Niemack et al. 2010) and SPTpol 
(McMahon et al. 2009) coming online in the near future. 

With continued SZ surveys such as ACT and SPT and 
their polarization counterparts, in addition to data forthcom- 
ing from the Planck satellite, we will no doubt increase the 
number of SZ cluster detections. We anticipate that upcom- 
ing larger galaxy cluster catalogs will make significant contri- 
butions to our understanding of both cluster astrophysics and 
cosmology. 
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APPENDIX 


A. SELF-SIMILAR SCALING RELATION BETWEEN SZ SIGNAL AND MASS 
If clusters were self-similar and isothermal, then we would expect the scaling relation between SZ signal and mass to be 

Lhaio « M^ o E(z) 2 ' n f^/4 , ( A 1 ) 

where Lhaio is the Compton y-parameter integrated over the surface of the cluster in units of aremin 2 , and E(z) = [Q m (l + z) 3 + 
Oa ] 1/2 for a flat ACDM cosmology. The angular diameter distance is denoted by d A , and ,/g as is the gas mass fraction. For the 
Compton y-parameter integrated over a fixed aperture we have 

WeOcLhaiof^^V, (A2) 

\ ^halo / 


where Raperture oc d A and Ahaio oc Ml!lE(z) 2/3 . Note that this equation is appropriate if the aperture size is smaller than the size 
of the cluster. The above gives 


/aperture X A/halo E(z) /gas- 


(A3) 


To write L aperturc as a function of (1 +z) we note that at z=0.5 (the mean redshift of our cluster sample) E(z) oc (1 +z) 0 - 835 for 
Q m = 0.27. Thus 


/aperture OC 47^ ajo ( 1 +z) , 


(A4) 


where B = 1 .0 and C= 1.67. 
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B. CLUSTER LIKELIHOOD FUNCTION 

Below we describe the construction of the likelihood function for SZ clusters detected in millimeter-wave surveys. From 
Poisson statistics, the probability of observing «, counts expecting A, counts is 


Print |A,) = 


A"-g- A - 

«j! 


(Bl) 


Given a data set of {«,} counts in Ni, observed bins and a corresponding prediction, {A/}, the probability of the data given the 
prediction is 


i=i 


A”'<T a ' 


»/! 


(B2) 


where A,- = P r (y** ,A,B,C,S, {c/})!VAy obs Az obs . Here P r (y ob \z obs ,A,B,C,S, {c y }) is the probability of observing a cluster in 
bin i, and A is a normalization factor giving A* units of counts (see below). The observed SZ signal and redshift of a given cluster 
are denoted by y obs and ~ obs , and Ay obs and Az obs denote the size of the bin. The parameters A,B,C,S describe the scaling relation 
between SZ signal and mass and are defined below. The cosmological parameters are indicated by { Cj } . 

If we allow the bin sizes to be small enough that each observed bin holds no more than one observed cluster, then 


n N/> 

lnPr({«/}|{A,}) = InA/ A,- = ln£ (B3) 

M j=l 

as given in Cash (1979). Note that the In (w,!) term has been dropped as it is independent of any change in parameters, and n 
represents the total number of clusters observed. Thus we have 


ln£ = Y ln(P r (y obs ; z obs .A,B,C,S, {cj})Ndy obs dz obs )- I I P r (y obs ,z ob \A,B.C,S , { C j})Ndy^di 


r obs 


with 


P r (y obs , z obs ,A.B,C,S. { C] }) = jjj P r (y obs , z obs ,A,B,C,S, {eyD/Az^nA^D^ 


P,.(y obs . z obs \A,B,C.S, {cy} ,/ rae . z 1 ™', In A/ m )P r (A ,B,C,S, { cy } ,y true ,z true , \nM ime )dy tmc dz ttue d\nM ln 


(B4) 


(B5) 

(B6) 


P r (y ob \z ob ^ e ,z^)P r (y tme \A,B,C,S,{cj},z^^\n^ me )P r (A,B,C,S,{cj} ! z^ e MAf Me )dy^ e dz^dlnAi tme (B7) 


/* oo pOO pOO 

= / d\nirf me / dy m,e / dz tme P r (y ob!< [y 1 rac )P r (z obs Iz 1 ™ 6 )P r (y tm<: \A.B,C,S, z trae , lnjl/ rae )P r (lnM true , z 1 ™' | { c ; } ) 

J- oo JO Jo 

xP r ({cj} )P r (A )P r (B)P r (C)P r (S) 

using the definition of conditional probability. Here P r ({c/}) is any external prior on {cy} such as a WMAP prior. 

We assume the following SZ signal - mass scaling relation with log normal scatter, S, 

,/r'\ ! /l +z bC 


A M o ) 


1 +2 


(B8) 


(B9) 


where Mq- 5 x 10 l4 M Q h 1 and zq = 0.5. This gives 


Pr(y lme \A,B,C,S, z ,rue , In A/ rae ) = 


V2nSy" 


■ exp 


-(lnv true - BlnA/"' llc - Cln( 1 +z true ) - lrr4 + BlnM 0 + Cln( 1 +z 0 )) 2 

2S 2 


We also assume Gaussian priors on the scaling relation parameters as indicated by simulations, giving 


P r (A) - 


1 


\Z2tt <t\ 


exp 


-(A-Ao ) 1 

2a\ 


Similar relations hold for P r (B), P r (C ), and P r (S). 
From the mass function we have 


(BIO) 


(Bll) 
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P r (lnM lrue ,z true |{cy}) = 


{cj}) dV(z lme , {c ; }) 1 
r/lnM me dz tme N’ 


(B12) 


where here n is the number density of clusters. N is the total number of clusters when the above mass function is integrated over 
c/lniW true and dz tme . 

We also assume for the uncertainty on the observed SZ signal and redshift that 


P r(y obs |/ uc )=- 
P,.(z obs |z ,nle ) = 


l 

V2ix (Jobs 

1 

\/2tto 


exp 


- exp 


_^y>bs_yrue)2 

2<J obs 


-(z o1 


ibs _true\2 




(B 1 3) 

(B14) 


where these two expressions should also be multiplied by l+erf(v i r ^- / /^ 1 g 0 . ] s ' nce the limits of integration are from 0 to oo. 
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